Degenerate atom-molecule mixture in a cold Fermi gas 
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We show that the atom-molecule mixture formed in a degenerate atomic Fermi gas with in- 
terspecies repulsion near a Feshbach resonance, constitutes a peculiar system where the atomic 
component is almost non-degenerate but quantum degeneracy of molecules is important. We de- 
velop a thermodynamic approach for studying this mixture, explain experimental observations and 
predict optimal conditions for achieving molecular BEC. 
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Interactions between particles play a crucial role in the 
behavior of degenerate quantum gases. For instance, the 
sign of the effective mean field interaction determines the 
stability of a large Bose-Einstein condensate (BEC), and 
the shape of such a condensate in a trap can be signifi- 
cantly altered from its ideal gas form y). In degenerate 
Fermi gases the effects of mean field interactions are usu- 
ally less pronounced in the size and shape of the trapped 
cloud, and these quantities are mostly determined by 
Fermi statistics. The strength of the interactions, how- 
ever, can be strongly increased by making use of a Fesh- 
bach resonance 0,11] j and then the situation changes. 

Recent experiments present two types of measurement 
of the interaction energy in a degenerate two-component 
Fermi e;as near a Feshbach resonance 

0, II H la At 

JILA [j| and at MIT the mean field energy was found 
from the frequency shift of an RF transition for one of 
the atomic states. The results are consistent with the 
magnetic field dependence of the scattering length a, the 
energy being positive for a > and negative for a < 0. 
In the Duke 3 and ENS 0] experiments with 6 Li, the 
results are quite different. The interaction energy was 
obtained from the measurement of the size of an expand- 
ing cloud released from the trap. A constant ratio of the 
interaction to Fermi energy, Ei nt /Ep rts —0.3, was found 
around resonance, irrespective of the sign of a 0,0- R 
was explained in Ref. by claiminga universal behavior 
in this strongly interacting regime yj. The ENS studies 
in a wide range of magnetic fields |6| found that E mt 
changes to a large positive value when a is tuned posi- 
tive, but only at a field strongly shifted from resonance. 

In contrast to the JILA and MIT Q studies provid- 
ing a direct measurement of the mean field interaction en- 
ergy, the Duke £| and ENS experiments measure the 
influence of the interactions on the gas pressure. An in- 
terpretation of the ENS experiment involves the creation 
of weakly bound molecules via three-body recombination 
at a positive a Q. Far from resonance, the binding en- 
ergy of the produced molecules and, hence, their kinetic 
energy are larger than the trap depth and the molecules 
escape from the trap. The interaction energy is then 
determined by the repulsive interaction between atoms 
and is positive Close to resonance, the three-body 



recombination is efficient |9j and the molecules remain 
trapped as their binding energy es is smaller than the 
trap depth 0, 0. They come to equilibrium with the 
atoms, reducing the pressure in the system. 

Away from resonances, the interaction strength is pro- 
portional to a, and is given by g = 4nK a/M, with M the 
atom mass. Close to resonance this relation is not valid, 
as the value of |o| diverges to infinity and the scatter- 
ing process strongly depends on the collision energy. For 
Boltzmann gases, already in the 1930's, Beth and Uh- 
lenbeck calculated the second virial coefficient by in- 
cluding both the scattering and bound states for the rel- 
ative motion of pairs of atoms [Tl| . A small interaction- 
induced change of the pressure in this approach is nega- 
tive on both sides of the resonance [fH [llj • 

However, current experiments are not in the Boltz- 
mann regime. In this letter we show that the atom- 
molecule mixture formed in a cold atomic Fermi gas, con- 
stitutes a peculiar system in which the atomic component 
is almost non-degenerate, whereas quantum degeneracy 
of the molecules can be very important. This behavior 
originates from a decrease of the atomic fraction with 
temperature. It is present even if the initial Fermi gas 
is strongly degenerate in which case almost all atoms are 
converted into molecules. We develop a thermodynamic 
approach for studying this mixture, predict optimal con- 
ditions for achieving molecular BEC, and properly de- 
scribe the interaction effects as observed at ENS p. 

We assume that fermionic atoms are in equilibrium 
with weakly bound (bosonic) molecules formed in the re- 
combination process. The molecules are treated as point 
bosons. Atom-molecule and molecule-molecule interac- 
tions are omitted at first, and will be discussed in a later 
stage. For a large scattering length a > 0, the binding 
energy of the weakly bound molecules is = ?i 2 / (Ma 2 ), 
and their size is roughly given by a/2. For treating them 
as point bosons, this size should be smaller than the 
mean interparticle separation. This requires the inequal- 
ity n(a/2) 3 < I, which at densities n ~ I0 13 cm -3 is 
satisfied for a < I8000ao, and excludes a narrow vicinity 
of the Feshbach resonance. 

The presence of molecules reduces the number of par- 
ticles in the atomic component and to an essential extent 
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FIG. 1: Fraction of unbound atoms n a /n (lower curves, bold) 
and fraction of atoms bound into molecules, 2n m /n, (upper 
curves) versus TVes: a) nAy = 2.5, squares and circles show 
the ENS data 12011 : b) nAy = 14.8, and the vertical lines indi- 
cate the onset of molecular BEC. Dashed curves are obtained 
including atom-molecule and molecule-molecule interactions. 



lifts its quantum degeneracy. The molecular chemical po- 
tential is negative in the absence of atom-molecule and 
molecule-molecule mean field, and thermal equilibrium 
between atoms and molecules requires a negative chemi- 
cal potential for the atoms. We thus assume a priori that 
the occupation numbers of the states of atoms are small. 
This proves to be the case at any temperature, except 
for very low T where the atomic fraction is negligible. 
Under these conditions we omit pairing correlations be- 
tween the atoms, which are important for describing a 
crossover from the BCS to BEC regime 0, 0, 0, 0] 
and can be expected even in the non-superfluid state. 

Assuming equal densities of the atomic components, 
labeled as j and j, their chemical potentials arc fi-\ — 
Hi = fi, where \i is the chemical potential of the system 
as a whole. The molecular chemical potential is fi m = 
— €b + firm with fl m < being the chemical potential of 
an ideal gas of bosons with the mass 2M. The condition 
of thermal equilibrium, /if + /if = /j, m , then reads 



2/j, 



(1) 



From Eq. we will obtain the number of molecules N m 
and the number of atoms N a for given temperature T and 
total number of atomic particles TV = N a + 2N m . This 
requires us to obtain the expression for the occupation 
numbers of the atoms and the dependence of /i on N a . 

The main difficulty with constructing a thermody- 
namic approach for the degenerate molecule-atom mix- 
ture is related to the resonance momentum-dependent 
character of the atom-atom interactions. This difficulty is 
circumvented for small occupation numbers of the atoms. 
Then, even at resonance, the interaction energy is equal 
to the mean value of the interaction potential for a given 



relative momentum of a colliding pair, averaged over the 
momentum distribution. In this respect, the interaction 
problem becomes similar to the calculation of the total 
energy of a heavy impurity as caused by its interactions 
with the surrounding electrons in a metal |l9l |. This ap- 
proach leads to a relation between the collision-induced 
shift of the energy levels of particles in a large spherical 
box, and the scattering phase shift. Adding the integra- 
tion over the states of the center of mass motion for pairs 
of atoms, we find that the total energy of interatomic in- 
teraction is equal to ^Zkk' Skk'^T^' A 4 ; T)vi(k' , fi, T)/V, 
where and are occupation numbers of sin gle-p article 
momentum states, and V is the volume (cf. [H|). The 
momentum-dependent coupling constant is given by 
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The phase shift 8 is expressed through the relative mo- 
mentum q = |k — k'l/2 and the scattering length a as 
S = — arctanqa. In the limit of q\a\ <C 1, Eq.® trans- 
forms into the ordinary coupling constant g = 4irh 2 a/M. 

As we have i/f(fc,/x,T) = v^ik, (j,,T) = v^, the total 
energy of the atomic component and the number of par- 
ticles in this component can be written in the form 
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N a = 2j2"k- (3) 



In our mean-field approach, the entropy of the atoms is 
given by the usual combinatorial expression |l8j : 

S a = -2^Kln^ fe + (1 - ^)ln(l - v k )]. (4) 



Equations J2J and 0} immediately lead to an expression 
for the atomic grand potential Q a = E a — TS a — [iN a . 
Then, using the relation A^ a = —(dfl a /d[i)T.v, we obtain 
for the occupation numbers of atoms: 



v k = [exp{(e fe - M )/T} + l]- 



(5) 



where e k = h k 2 /2M + U k , and U k = J2k' 9wVk>/V is 
the mean field acting on the atom with momentum k. 
Accordingly, the expression for the grand potential and 
pressure of the atomic component reads: 



tt a = -P a V = y>Tln(l - v k ) - Ukv k \- 
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This set of equations is completed by the relation be- 
tween the density of bosonic molecules and their chemical 
potential. In the absence of molecular BEC we have: 



(V2/A T ) 3/2 ff 3 /2(cxp(A m /T), 



(7) 



where g a (x) = x j /j a , and A T = (2nh z /MT) 1 / 2 is 

the thermal de Broglie wavelength for the atoms. For 
n m Ay > 7.38 the molecular fraction becomes Bosc- 
condensed, and we have jl = and fi = —eb/2. Similarly, 
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FIG. 2: Molecular degeneracy parameter n m A^ under adia- 
batic variation of a for 6 Li, assuming nA^ = 15 close to reso- 
nance. The dashed curve is obtained including atom-molecule 
and molecule-molecule interactions. The horizontal dashed 
line shows the critical value for molecular BEC. 



the energy, entropy, and grand potential of the molecules 
are given by usual equations for an ideal Bose gas [TlT ] . 

From Eqs. CEJ-Q) we obtain the fraction of unbound 
atoms n a /n and the fraction of atoms bound into 
molecules, 2n m /n, as universal functions of two param- 
eters: T/eg and nA T , where n is the total density of 
atomic particles. The dependence of atomic and molec- 
ular fractions on T/eg for two values of nA T is shown in 
Figffl The molecular fraction increases and the atomic 
fraction decreases with decreasing T/eg. Occupation 
numbers of the atoms are always small, whereas quan- 
tum degeneracy of molecules is important. The dotted 
line in Fig. lb indicates the onset of molecular BEC. 

This mixture was realized in the ENS experiment 0, 
where the occupation numbers for the molecules were 
up to 0.3 and the molecular fraction was exceeding the 
atomic one. In the recent studies [53, 0, H3, H3 al- 
most all atoms were converted into molecules by swee ping 
the magnetic field across the resonance, and at ENS [2£| 
the temperature was within a factor of 2 from molecular 
BEC. Remarkably, one can modify the molecular fraction 
and degeneracy parameter n m A T by adiabatically tuning 
the atom-atom scattering length, as shown in Fig. 2. The 
decrease of a increases the binding energy cr and the 
molecular fraction, and thus causes heating [2(|. Close 
to resonance, n m A T remains almost constant and then 
decreases due to heating. 

The atom-molecule and molecule-molecule interactions 
are readily included in our approach for a <C At, where 
the corresponding coupling constants are g am = 0.9g 
and g m = 0.3g pr3 |. In this limit the interactions 
provide an equal shift of the chemical potential and 
single-particle energy e^. For the atoms this shift is 
n a g/2 + n m g am , where the first term is the atom-atom 
contribution [/&. For the (non-condensed) molecules the 
shift is n a g am + 2n m g m . The entropy of the mixture is 
given by the same expressions as in the absence of the in- 
teractions. As seen in Fig.l and Fig. 2, the atom-molecule 
and molecule-molecule interactions do not significantly 
modify our results. From Fig. 2 one then concludes that 
the conditions for achieving molecular BEC arc optimal 



for values of a as low as possible while still staying at 
the plateau, as at larger a the interaction between the 
molecules can reduce the BEC transition temperature 0] . 

We now analyze the interaction effect observed at ENS 
for trapped clouds in the hydrodynamic regime The 
experiment was done near the Feshbach resonance lo- 
cated at the magnetic field Bq = 810 G, and the data 
results from two types of measurements of the size of the 
cloud released from the optical trap. In the first one, 
the magnetic field and, hence, the scattering length, are 
kept the same as in the trap. Therefore, the cloud ex- 
pands with the speed of sound c s = (dP/dp)s, where 
p = ran is the mass density. The speed c s and, hence, the 
size of the expanding cloud are influenced by the presence 
of molecules and by the interparticle interactions. 

In the second type of measurement, the magnetic field 
is first rapidly ramped down and the scattering length 
becomes almost zero on a time scale t ~ 2as. This time 
scale is short compared to the collisional time. Therefore, 
the spatial distribution of the atoms remains the same as 
in the initial cloud, although the mean field is no longer 
present. At the same time, a rapid decrease of a increases 
the binding energy of molecules eg. However, as the time 
t < h/eg, they can not adiabatically follow to a deeper 
bound state and dissociate into atoms which acquire ki- 
netic energy. Thus the system expands symmetrically as 
an ideal gas of N atoms, with the initial density pro- 
file. The momentum distribution fk will be a sum of 
the initial atomic momentum distribution and one that 
arises from the dissociated molecules. The latter is found 
assuming an abrupt change of o and, hence, projecting 
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FIG. 3: Calculated (solid line) and measured (squares 
and crosses) ratio /3 of the interaction to kinetic energy (see 
text). The calculated line for B > 790G is for experimental 
conditions T = 0.9E F = 3.4/iK and n = 3 x 10 13 cm~ 3 . For 
B < 700G we take the averaged experimental conditions T = 
1.1E F = 2.4^iK and n = 1.3 x 10 13 cm -3 . For 700 < B < 
800G, we use the local conditions (see 0]). Inset: Scattering 
length as a function of magnetic field. 
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the molecular wave function on a complete set of plane 
waves. This gives rise to a distribution c{q) for the rela- 
tive momentum q. The single-particle momentum distri- 
bution for the atoms produced out of molecules results 
then from convoluting c(|k — k'|/2) with the molecular 
distribution function v m (k + k') by integrating over k'. 
One can establish a relation between the expansion veloc- 
ity Vq of this non-equilibrium system and the expansion 
velocity cq of an ideal equilibrium two-component atomic 
Fermi gas which has the same density and temperature: 
47r 3 n fM-oo/h fikk 2 f k = J^ 00 dkk 2 v k , with v k being the 
ideal-gas momentum distribution. Using the scaling ap- 
proach 0, HH , one can find that in the spherical case 
the velocity Co coincides with the expansion velocity of 
the hydrodynamic Fermi gas in the absence of mean-field 
interactions and, accordingly, is given by Cq = hPo/Zp, 
where Pq = 2E /3V is the pressure. 

The relative difference between the squared size of the 
expanding cloud in the two described cases can be treated 
as the ratio of the interaction to kinetic energy and called 
the interaction shift. This interaction shift is then given 
by the relative difference between the two squared veloc- 
ities: f3 = [c 2 — Vq]/v$. Our results for this quantity are 
calculated for experimental conditions and are presented 
in Fig|3 The sound velocity c s was obtained using the 
above developed approach including only atom-atom in- 
teractions. The field region where n(a/2) 3 > 1, is beyond 
the validity of this approach and is shown by the dashed 
curve. In Fig|3|we also show our previous results for fields 
B > 810 G (a < 0) and B < 700 G (0 < a < 2000a ), 
where molecules are absent @. 

Our quantum-statistical approach gives a negative in- 
teraction shift on both sides of the Feshbach resonance, in 
good quantitative agreement with the experiment. With- 
out molecules present, the interaction energy would jump 



to positive values left from resonance, as can be seen from 
our calculation in Ref. @. This demonstrates that the 
apparent field shift from resonance, where a sign-change 
in the interaction energy is observed, is an indirect sig- 
nature of the presence of molecules in the trap. 

For high temperatures T 3> Ep and small binding en- 
ergy eb *C T, we find that (3 has a universal behavior 
and is proportional to the second virial coefficient. How- 
ever, this only holds at high temperatures (cf. 01), and 
at low T the molecule-molecule interaction can strongly 
influence the result. For T approaching the temperature 
of molecular BEC, which is T c fx h 2 n 2 / 3 /M w 0.2E F , the 
atomic fraction is already small and the sound velocity 
c s is determined by the molecular cloud. For a -C At we 
find c 2 = 0AT c /M + ng m /2M, where the second term 
is provided by the molecule-molecule interaction and is 
omitted in the high-T approach. The ratio of this term 
to the first one is ~ 5(na 3 ) 1 / 3 . For B = 700 G at den- 
sities of Ref. , it is equal to 1 and is expected to grow 
when approaching the resonance. 

Thus, except for a narrow region where n|a| 3 3> 1, 
one can not speak of a universal behavior of the shift f3 
on both sides of the resonance. The situation depends 
on possibilities of creating an equilibrium atom-molcculc 
mixture. Moreover, at low temperatures the universality 
can be broken by the molecule-molecule interactions. 
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